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Hard Constraints and the Bethe Lattice: 
Adventures at the Interface of 
Combinatorics and Statistical Physics 

Graham R. Brightwell* Peter Winkler^ 

Abstract 

Statistical physics models with hard constraints, such as the discrete hard- 
core gas model (random independent sets in a graph), are inherently combi- 
natorial and present the discrete mathematician with a relatively comfortable 
setting for the study of phase transition. 

In this paper we survey recent work (concentrating on joint work of the au- 
thors) in which hard-constraint systems are modeled by the space Hom(G, H) 
of homomorphisms from an infinite graph G to a fixed finite constraint graph 
H. These spaces become sufficiently tractable when G is a regular tree (of- 
ten called a Cayley tree or Bethe lattice) to permit characterization of the 
constraint graphs H which admit multiple invariant Gibbs measures. 

Applications to a physics problem (multiple critical points for symmetry- 
breaking) and a combinatorics problem (random coloring), as well as some 
new combinatorial notions, will be presented. 
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1. Introduction 

Recent years have seen an explosion of activity at the interface of graph theory 
and statistical physics, with probabilistic combinatorics and the theory of computing 
as major catalysts. The concept of "phase transition", which a short time ago 
most graph theorists would barely recognize, has now appeared and reappeared in 
journals as far from physics as the Journal of Combinatorial Theory (Series B). 
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Traffic between graph theory and statistical physics is already heavy enough 
to make a complete survey a book-length proposition, even if one were to assume a 
readership with knowledge of both fields. 

This article is intended for a general mathematical audience, not necessarily 
acquainted with statistical physics, but it is not to serve as an introduction to the 
field. Readers are referred to texts such as El Q3| f° r more background. 

We will present only a small part (but we hope an interesting one) of the interface 
between combinatorics and statistical physics, with just enough background in each 
to make sense of the text. We will focus on the most combinatorial of physical 
models — those with hard constraints — and inevitably on the authors' own research 
and related work. 

We hope it will be clear from our development that there is an enormous 
amount of fascinating mathematics to be uncovered by studying statistical physics, 
quite a lot of which has been or will be connected to graph theory. What follows is 
only a sample. 

2. Random independent sets 

In what follows a graph G = (V, E) consists of a set V (finite or countably 
infinite) of nodes together with a set E of edges, each of which is an unordered pair 
of nodes. We will sometimes permit loops (edges of the form {v,v}) but multiple 
edges will not be considered or needed. We write u ~ v, and say that u is "adjacent" 
to v, if {u, v} S E; a set U C V is said to be independent if it contains no edges. 

The degree of a node u of G is the number of nodes adjacent to u; all graphs 
considered here will be locally finite, meaning that all nodes have finite degree. A 
path in G (of length k) is a sequence Uo, u\, . . . , Uk of distinct nodes with m ~ Ui+i] 
if in addition Uk ~ uq we have a cycle of length fe+1. If every two nodes of G are 
connected by a path, G is said to be connected. 

The plane grid Z 2 is given a graph structure by putting (i,j) ~ iff 
K — *'l + = 1- Let us carve out a big piece of Z 2 , say the box f? 2 := 

{{hj) I — 71 — h i — n }- Let 7 be a uniformly random independent set in the graph 
-B 2 ; in other words, of all sets of nodes (including the empty set) not containing an 
edge, choose one uniformly at random. What does it look like? 

Plate 1 shows such an I (in a rectangular region). Here the nodes are repre- 
sented by squares, two being adjacent if they have a common vertical or horizontal 
border segment. The sites belonging to I are colored, the others omitted. It is by 
no means obvious how to obtain such a random independent set in practice; one 
cannot simply choose points one at a time subject to the independence constraint. 
In fact the set in Plate 1 was generated by Markov chain mixing, an important and 
fascinating method in the theory of computing, which has by itself motivated much 
recent work at the physics-combinatorics interface. 

It is common in statistical physics to call the nodes of _B 2 sites (and the edges, 
bonds). Sites in I are said to be occupied; one may imagine that each occupied site 
contains a molecule of some gas, any two of which must be at distance greater than 
1. 
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In the figure, even occupied sites (nodes € I for which i+j = mod 2) 
are indicated by one color, odd sites by another. A certain tendency for colors to 
clump may be observed; understandably, since occupied sites of the same parity 
may be as close as \/2 (in the Euclidean norm) but opposite-parity particles must 
be at least apart. 

It stands to reason that if more "particles" were forced into /, then we might 
see more clumping. Let us weight the independent sets according to size, as follows: 
a positive real A, called the activity (or sometimes fugacity) is fixed, and then each 
independent set / is chosen with probability proportional to \W. We call this the 
"A-measure" . Of course if A = 1 we are back to the uniform measure, but if A > 1 
then larger independent sets are favored. 

The A-measure for A ^ 1 is, to a physicist, no less natural than the uniform. 
In a combinatorial setting, such a measure might arise e.g. if the particles happen 
to be of two different types, with all "typed" independent sets equiprobable; then 
the probability that a particular set of sites is occupied is given by the A-measure 
with A = 2. 

Plate 2 shows a random / chosen when A = 3.787. The clusters have grown 
hugely as more particles were packed in. Push A up just a bit more, to 3.792, and 
something like Plate 3 is the result: one color (parity) has taken over, leaving only 
occasional islands of the other. 

Something qualitative has changed here, but what exactly? The random in- 
dependent sets we have been looking at constitute what the physicists call the 
hard-core lattice gas model, or "hard-core model" for short. Readers are referred to 
the exceptionally readable article |2j in which many nice results are obtained for this 
model 1 . On the plane grid, the hard-core model has a "critical point" at activity 
about 3.79, above which the model is said to have experienced a phase transition. 

3. What is phase transition? 

There is no uniformity even among statistical physicists regarding the defi- 
nition of phase transition; in fact, there is even disagreement about whether the 
"phases" above are the even-dominated versus odd-dominated configurations at 
high A, or the high-A regime versus the low. Technical definitions involving points 
of non-analyticity of some function miss the point for us. 

The point really is that a slight change in a parameter governing the local 
behavior of some statistical system, like the hard-core model, can produce a global 
change in the system, which may be evidenced in many ways. For example, suppose 
we sampled many independent sets in B\ at some fixed A, and for each computed 
the ratio of the number of even occupied sites to the number of odd. For low A 
these numbers would cluster around 1/2, but for high A they would follow a bimodal 
distribution; and the larger the box size n, the sharper the transition. 

Here's another, more general, consideration. Suppose we look only at inde- 
pendent sets which contain all the even sites on the boundary of B\. For these / 

1 Readers, however, are cautioned regarding conducting a web search with key-words "hard- 
core" and "model". 



608 



G. R. Brightwcll P. Winkler 



the origin would be more likely to be occupied than, say, one of its odd neighbors. 
As n grows, this "boundary influence" will fade — provided A is low. But when A is 
above the critical point, the boundary values tend to make / an even-dominated set, 
giving any even site, no matter how far away from the boundary, a non-disappearing 
advantage over any odd one. 

Computationally-minded readers might be interested in a third approach. Sup- 
pose we start with a fixed independent set Iq, namely the set of all even sites in B^, 
and change it one site at a time as follows: at each tick of a clock we choose a site it 
at random. If any of it's neighbors is occupied, we do nothing. Otherwise we flip a 
biased coin and with probability A/(l + A) we put it in / (where it may already have 
been), and with probability 1/(1 + A) we remove it (or leave it out). The result is a 
Markov chain whose states are independent sets and whose stationary distribution, 
one can easily verify, is exactly our A-measure. Thus if we do this for many steps, 
we will have a nearly perfect sample from this distribution — but how many steps 
will that take? We believe that when A is below its critical value, only polynomi- 
ally (in n) steps are required — the Markov chain is said to be rapidly mixing] even 
polylogarithmic, if we count the number of steps per site. But for high A it appears 
to take time exponential in n (or perhaps in ^/n) before we can expect to see an 
odd-dominated independent set. The exact relationship between phase transition 
and Markov chain mixing is complex and the subject of much study. 

All these measures rely on taking limits as the finite box _B 2 grows; the very 
nice discovery of Dobrushin, Lanford and Ruelle |U E3 is that there is a way to un- 
derstand the phenomenon of phase transition as a property of the infinite plane grid. 
The idea is to extend the A-measure to a probability distribution on independent 
sets on the whole grid, then ask whether the extension is unique. 

We cannot extend the definition of the A-measure directly since A' 7 ' is generally 
infinite, but we can ask that it behave locally like the finite measure. We say that a 
probability distribution /i on independent sets in the plane grid is a Gibbs measure 
if for any site it the probability that it is in /, given the sites in I D (Z 2 \ {it}), is 
A/(l + A) if the neighborhood of it is unoccupied and, of course, otherwise. 

It turns out that Gibbs measures always exist (here, and in far greater gener- 
ality) but may or may not be unique. When there is more than one Gibbs measure 
we will say that there is a phase transition. For the hard-core model on Z 2 , there 
is a unique Gibbs measure for low A; but above the critical value, there is a Gibbs 
measure in which the even occupied sites are dominant and another in which the 
odd sites are dominant (all other Gibbs measures are convex combinations of these 
two). How can you construct these measures? Well, for example, the even measure 
can be obtained as a limit of A-measures on boxes whose even boundary sites are 
forced to be in /. The fact that the boundary influence does not fade (in the high 
A case) implies that the even and odd Gibbs measures are different. 

We have noted that the critical value of A for the hard-core model on 1? is 
around 3.79. This is an empirical result and all we mathematicians can prove is 
that there is at least one critical point, and all such are between 1.1 and some high 
number. It is believed that, for each d, there is just one critical value A^ on Z d . 
It is also to be expected that A^ is decreasing in d, but only recently has it been 
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shown that the largest critical value on Z d tends to as d — > oo. This result was 
obtained by David Galvin and Jeff Kahn two combinatorialists, using graph 
theory, geometry, topology, and lots of probabilistic combinatorics. A consequence 
of their work is that A = 1 is above the critical value(s) for sufficiently large d; this 
can be stated in a purely combinatorial way: for sufficiently high d and large n, 
most independent sets in are dominated by vertices of one parity. 

In the next section we explain how we can use graphs to understand models 
with hard constraints; then, in the section following that, we will switch from Z 2 to 
a much easier setting, in which we can get our hands on nice Gibbs measures. 

4. Hard constraints and graph homomorphisms 

We are interested in what are sometimes called "nearest neighbor" hard con- 
straint models, where the constraints apply only to adjacent sites. Each site is to be 
assigned a "spin" from some finite set, and only certain pairs of spins are permitted 
on adjacent sites. We can code up the constraints as a finite graph H whose nodes 
are the spins, and whose edges correspond to spins allowed to appear at neighboring 
sites. This constraint graph H may have some loops; a loop at node v € H would 
mean that neighboring sites may both be assigned spin v. We adopt the statistical 
physics tradition of reserving the letter u q" for the number of spins, that is, the 
number of nodes in H. 

The graph G (e.g. Z 2 , above) of sites, usually infinite but always countable 
and locally finite, is called (by us) the board. A legal assignment of spins to the 
sites of G is nothing more or less than a graph homomorphism from G to H, i.e. a 
map from the sites of G to the nodes of H which preserves edges. We denote the 
set of homomorphisms from G to H by Hom(G, H) , and give it a graph structure 
by putting ip ~ tp if ip and ip differ at exactly one site of G. 

We will often confuse a graph with its set of nodes (or sites). In particular, if 
U is a subset of the nodes of G then U together with the edges of G contained in 
U constitute the "subgraph of G induced by U" , which we also denote by U. 

In the hard-core model, the constraint graph H consists of two adjacent nodes, 
one of which is looped: a function from a board G to this H is a homomorphism iff 
the set of sites mapped to the unlooped node is an independent set. Plate 4 shows 
some constraint graphs found in the literature. 

When H is complete and every node is looped as well, there is no constraint 
and nothing interesting happens. 

When H is the complete graph K q (without loops), homomorphisms to H are 
just ordinary, "proper" g-colorings of the board. (A proper q-coloring of a graph 
G is a mapping from the nodes of G to a g-element set in which adjacent nodes 
are never mapped to the same element.) This corresponds to something called the 
"anti-ferromagnetic Potts model at zero temperature" . In the anti-ferromagnetic 
Potts model at positive temperature, adjacent sites are merely discouraged (by an 
energy penalty), not forbidden, from having the same spin; thus this is not a hard 
constraint model in our terminology. The q — 2 case of the Potts model is the 
famous Ising model. 
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For a general constraint graph H, we need to elevate the notion of activity to 
vector status. To each node i of H we assign a positive real activity A$, so that 
H now gets an activity vector A := (Ai, . . . , X q ). When the board G is finite, each 
homomorphism ip € Hom(G, H) is assigned probability proportional to 

n ■ 

v£G 

We can think of A^ as the degree to which we try to use spin i, when it is available. 
For example, if we know the spins of the neighbors of site v and consequently, 
say, spins i, j and k are allowed for v, then the A-measure forces Pr(<^(t>) = i) = 
Xi/iXi + Xj + Xk)- 

When G is infinite, things get a little more complicated. A finite subset (and 
its induced subgraph) U C G will be called a "patch" and its boundary dU is the 
set of sites not in U but adjacent to some site of U. We define U + := U U dU. If (p 
is a function on G, then ip \ U denotes its restriction to the subset U . 

We say that /i is a Gibbs measure for A if: for any patch U C G, and almost 
every ip e Hom(G,i/), 

Pr(ip\U = ip\U\ip\(G-U) = ip \(G -U)) = Pr (<p\U = ip\U\tp\dU = tp\dU) 
fj. u+ 

where "Prj/+" refers to the finite A-measure on U + . 

This definition looks messy but it just means that the probability distribution 
of a random ip inside a patch U depends only on its value on the boundary of U, 
and is the same as if U and its boundary comprised all of the board. We will see 
later that when H has a certain nice property, as it does in the case of the hard-core 
model, it suffices to check the Gibbs condition only on patches consisting of a single 
site — we call this the one-site condition. 

It is a special case of a theorem of Dobrushin Jjjjj that there is always at least 
one Gibbs measure for any A on Hom(G, H); we are concerned with questions about 
when there is a unique Gibbs measure, and when there is a phase transition (i.e. 
more than one Gibbs measure). 

Let us again look briefly at possible implications for phase transition in the 
setting of finite boards. Given a finite board G, a constraint graph H and activities 
A, we define the point process V{G,H,X) as follows: starting from any element of 
Hom(G, H ), choose a site u of G uniformly at random, and give it a fresh spin ac- 
cording to the Gibbs condition, so that each 'legal' spin j is chosen with probability 
proportional to Xj. The point process is a Markov chain on Hom(G, H), and it is 
easy to check that the A-measure is a stationary distribution (which will be unique 
provided Hom(G, H) is connected, a point we will return to later). 

Running the point process for sufficiently long will thus generate a random 
homomorphism according to the A-measure. However, suppose that the finite board 
G is a large piece of an infinite board G' exhibiting a phase transition for our 
A. Then, if we start with a homomorphism arising from one Gibbs measure on 
Hom(G', H) (restricted to G), it is reasonable to expect that the point process will 
take a long time to reach a configuration resembling that from any other Gibbs 
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measure on Hom(G',i/). Thus it is generally believed that, in some necessarily 
loose sense, phase transition on an infinite graph corresponds to slow convergence 
for the point process on finite subgraphs. 

5. Cay ley trees and branching random walks 

Gibbs measures can be elusive and indeed it is generally a difficult task to prove 
that phase transitions occur on a typical board of interest, like Z d . In order to get 
results and intuition physicists sometimes turn to a more tractable board, called 
by them the Bethe lattice (after Hans Bethe) and by combinatorialists, usually, the 
Cayley tree. 

We denote by T r the r-branching Cayley tree, equivalcntly the unique con- 
nected (infinite) graph which is cycle- free and in which every site has degree r+1. 
T r is a vastly different animal from 7L d . It is barely connected, falling apart with 
the removal of any site; its patches have huge boundaries, comparable in size with 
the patch itself; its automorphism group is enormous. It's surprising that we can 
learn anything at all about Hom(Z d , H) from Hom(T r , H), and indeed we must be 
careful about drawing even tentative conclusions in either direction. Basic physical 
parameters like entropy become dodgy on non-amenable (big-boundary) boards like 
T r and a number of familiar statistical physics techniques become useless. More 
than making up for these losses, though, are the combinatorial techniques we can 
use to study Hom(T r , H). There are even situations (e.g. in the study of information 
dissemination) where T r is the natural setting. 

We are particularly interested in Gibbs measures on Hom(T r ,_ff) which have 
the additional properties of being simple and invariant. 

For any site u in a tree T, let d(u) be the number of edges incident with u and 
let Ci(u), C^u), . . . , Cm u ){u) be the connected components of T \ {u}. 

Definition 5.1 A Gibbs measure \x on Hom(T, H) is simple if, for any site u<ET 
and any node i € H, the ^-distributions of 

<P\Ci(u),...,<p\C d ( u) (u) 

are mutually independent given tp(u) = i. 

This condition, which is trivially satisfied by the A-measure for finite T, would 
follow from the Gibbs condition itself if fewer than two of the Ci(u)'s were infinite. 

Definition 5.2 Let A(G) be the automorphism group of the board G, and for any 
subset S C Hom(G, H) and k € A(G) let S o n := {<p o k : ^ £ S}. We say that 
a measure /i on Hom(G, H) is invariant if for any [i-measurable S C Hom(G, H) 
and any k £ A(G), we have fj,(S o k) = [i(S). 

Again, this condition is trivially satisfied for finite G; but for an infinite board 
with as many automorphisms as T r , it is quite strong. Later we consider relaxing 
it slightly. For now, we might well ask, how can we get our hands on any Gibbs 
measure for Hom(T r , H), let alone a simple, invariant one? 
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The absence of cycles in T r makes it plausible that we can get ourselves a Gibbs 
measure by building random configurations in Hom(T r ,iJ) one site at a time. We 
could choose a root x E T r , assign it a random spin i S H, then assign the neighbors 
of i randomly to the r+1 children of x; thereafter, each time a site u gets spin j we 
give its r children random spins from among the neighbors of j. 

The process we have described can be thought of as a branching random walk 
on H. Imagine amoebas staggering from node to adjacent node of H; each time an 
amoeba steps it divides into r baby amoebas which then move independently at the 
next time step. Of course, the (usually tiny) constraint graph H is shortly piled 
high with exponentially many amoebas, but being transparent they happily ignore 
one another and go on stepping and dividing. 

To get started we have to throw the first amoeba onto H where we imagine 
that its impact will cause it to divide r + 1 ways instead of the usual r. 

Note that the r = 1 case is just ordinary random walk, started somewhere on 
the doubly-infinite path T 1 and run both forward and backward. 

To determine what probabilities are used in stepping from one node of H to 
an adjacent node, we assign a positive real weight Wi to each node. For convenience 
we denote by Zi the sum of the weights of the neighbors of i (including i itself, if 
there is a loop at i). An amoeba-child born on node i then steps to node j with 
probability Wj/zi. If there is a loop at i, the amoeba stays at i with the appropriate 
probability, Wi/zi. 

Assuming H is connected and not bipartite 2 , the random walk (branching or 
not) will have a stationary distribution 7r; it is easily verified that 7Tj is proportional 
to WiZi for each i, and somewhat less easily verified that the mapping w c — > n is one- 
to-one provided w i nas been normalized to 1. We use the stationary distribution 
to pick the starting point for the first amoeba, i.e. to assign a spin to the root of 
T r . 

Finally, the payoff: not only does this node-weighted branching random walk 
give us a simple invariant Gibbs measure; it's the only way to get one. The following 
theorem appears in 0] but it is not fundamentally different from characterizations 
which can be found in Georgii |12| and elsewhere. 

Theorem 5.3 Let H be a fixed connected constraint graph with node-weights w 
and let r be a positive integer. Then the measure /i induced on Hom(T r ,i?) by the 
r-branching w-random walk on H is a simple, invariant Gibbs measure, for some 
activity A on H . Conversely, if H, r and X are given, then every simple, invariant 
Gibbs measure on Hom(T r ,i?) is given by the r-branching random walk on H with 
nodes weighted by some w. 

The proof is actually quite straightforward, and worth including here. Invari- 
ance of fi with respect to root-preserving automorphisms of T r is trivial, since the 
random walk treats all children equally; the only issue is whether the selection of 
root makes a difference. For this we need only check that for two neighboring sites 
u and v of T r , /1 is the same whether u is chosen as root or v is. But, either way 

2 A graph is bipartite if its nodes can be partitioned into two sets neither of which contains an 
edge. Thus, for example, the existence of a looped node already prevents H from being bipartite. 
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we may choose <p(u) and f(v) as the first two spins and the rest of the procedure 
is the same; so it suffices to check that for any (adjacent) nodes i and j of H, the 
probability that ip(u) = i and ip(v) — j is the same with either root choice. But 
these two probabilities are 

Wj Wi 

TTiPij = ZiWi = WiWj = ZjWj = ITjPji 

Zi Zj 

as desired. 

To show that \x is simple is, indeed, simple: if we condition on ip(u) = i then, 
using invariance to put the root at u, the independence of ip on the r+1 components 
of T r \ {u} is evident from the definition of the branching random walk. 

The activity vector A for which [i is a Gibbs measure turns out to be given by 

. _ w l 
K — —p ■ 

z 

Let U be any finite set of sites in T r , with exterior boundary dU. On account of 
invariance of labeling, we may assume that the root x does not lie in U + = J7U dU . 

Let g € Hom(?7 + , H); we want to show that the probability that a branching 
random walk ip matches g on U, given that it matches on dU, is the same as the 
corresponding conditional probability for the A-measure. 

Let T be the subtree of T r induced by U + and the root x; for any / e 
Hom(T, H), 

Pr(<p \T = f) = TT f{x) ■ ]J Pf(u)J(v) 



= Zf(x)w f{x ) n 



Z f(u) 



where u — > v means that v is a child of u in the tree. The factors Zf( u ) corresponding 
to sites u in U each occur as denominator r times in the above expression, since each 
site in U has all of its r successors in T; and of course each occurs once as a 

numerator as well. It follows that if we compare Pr(<£ \T = f) with Pr(</? \T = /'), 
where /' differs from / only on U, then the value of the first is proportional to 

n W f(u) _ TT \ 
-f— 11 X f(u) 

ueu /(«) ueu 

which means that /i coincides with the finite measure, as desired. 

Now let us assume that \i is a simple, invariant Gibbs measure on Hom(T r , H) 
with activity vector A, with the intent of showing that fi arises from a node- weighted 
branching random walk on H. 

We start by constructing a ^-random ip, site by site. Choose a root x of T 
and pick ip(x) from the a priori distribution a of spins of x (and therefore, by 
invariance, of any other site). We next choose a spin for the child y of x according 
to the conditional distribution matrix P = {pij} given by 

Pij := Pr (ip(y) = j \ tp(x) = i) ; 
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again, by invariance of /i, P is the same for any pair of neighboring sites. It follows 
that a — a ■ P, and moreover that P is the transition matrix of a reversible Markov 
chain, since the roles of x and y can be interchanged. 

Next we proceed to the rest of the children of x, then to the grandchildren, 
etc., choosing each spin conditionally according to all sites so far decided. 

We claim, however, that the distribution of possible spins of the non-root v 
depends only on the spin of its parent u; this is so because \i is simple and all sites 
so far "spun" are in components of T r \ {u} other than the component containing v. 
Thus the value of f(v) is given by P for every site d/i, and it follows that fj, arises 
from an r-branching Markov chain with state-space H , starting at distribution a. 

Evidently for any (not necessarily distinct) nodes i, j of H, there will be pairs 
(u, v) of adjacent sites with ip(u) — i and <p(v) — j if and only if i ~ j in H. Hence 
P allows transitions only along edges of H, and there is a unique distribution n 
satisfying n ■ P = 7r; thus a = tt. 

It remains only to show that P is a node- weighted random walk, and it turns 
out that a special case of the Gibbs condition for one-site patches suffices. Let j 
and j 1 be nodes of H which have a common neighbor i, and suppose that all of the 
neighbors of the root x have spin i. Such a configuration will occur with positive 
probability and according to the Gibbs condition for U = {x}, 

Pr(y(aQ = f) = Xf 
Pr&(x) = j) X, 

but 

Pr(</?(x) = i 

and similarly for j', so 



r+l 



Pr(^(a!) = j) n^t 1 



Thus the ratio 



Pj'i 



3 



is independent of i. 

Since P is reversible we have pij = itjPji , hence 

Pij/ ftj'Pj'i 



Pij TTjPji 

is also independent of i, and it follows that P is a node- weighted random walk on 
H . This concludes the proof of Theorem 15.31 

In view of Theorem 15. 31 if we can understand the behavior of the map w c — » A, 
we will know, given A, whether there is a nice Gibbs measure and if so whether 
there is more than one. The first issue is settled nicely in the following theorem, a 
proof of which can be found in 0] and requires some topology. A similar result was 
proved by Zachary |2.i|. 
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Theorem 5.4 For every r > 2, every constraint graph H and every set X of activ- 
ities for H , there is a node-weighted branching random walk on T" r which induces a 
simple, invariant Gibbs measure on Hom(T r , H). 

It's nice to know that we haven't required so much of our measures that they 
can fail to exist. 

A statistical physics dictum (true in great, but not unlimited, generality) says 
that there's never a phase transition in dimension 1; that holds here: 

Theorem 5.5 For any connected constraint graph H and any activity vector X, 
there is a unique simple invariant Gibbs measure on Hom(T 1 ,iJ). 

Furthermore, in any dimension, there's always some region where the map 
w A is one-to-one: 

Theorem 5.6 For any r and H there is an activity vector X for which there is only 
one simple invariant Gibbs measure on Hom(T r , H) . 

6. Fertile and sterile graphs 

The fascination begins when we hit an H and a A which boast multiple simple, 
invariant Gibbs measures. Let us examine a particular case, involving a constraint 
graph we call the "hinge" . 

The hinge has three nodes, which we associate with the colors green, yellow 
and red; all three nodes are looped and edges connect green with yellow, and yellow 
with red. Thus the only missing edge is green-red, and a (p € T r may be thought 
of as a green-yellow-red coloring of the tree in which no green site is adjacent to a 
red one. 

The hinge constraint in fact corresponds to a discrete version of the Widom- 
Rowlinson model, in which two gases (whose particles are represented by red and 
green) compete for space and are not permitted to occupy adjacent sites; see e.g. 
When A re d and A grce n are equal and large relative to A yo iiow, the Widom- 
Rowlinson model tends to undergo a phase transition as one gas spontaneously 
dominates the other. Plate 5 shows a red-dominated sample from the Widom- 
Rowlinson model on Z 2 , with the unoccupied sites left uncolored instead of being 
colored yellow. 

We can see the phase transition operate on the Cayley tree T 2 . If the green, 
yellow and red nodes are weighted 4, 2 and 1 respectively, A (normalized to integers) 
turns out to be (49, 18,49) — equal activity for green and red. How can a random 
walk which is biased so strongly toward green end up coloring a tree according 
to a Gibbs measure with symmetric specification? As a clue, let us examine a 
site u of T 2 which happens to be surrounded by yellow neighbors. To be colored 
green requires that a certain amoeba stepped from yellow to green, then both of 
its children returned to yellow. Thus the conditional probability that u is green is 
proportional to 
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as opposed to 



1 



( 



2 



) 



2 



4 + 2 + 1 



2 + 1 



for red, but these values are equal. 

Clearly the reversed weights 1, 2 and 4 would yield the same activity vector, 
and in fact a third, symmetric weighting, approximately 6, 7 and 6, does as well. 
Plate 6 shows pieces of T 2 colored according to these three weightings. Of course 
the colorings have different proportions and are easily identifiable; checking the 
stationary distributions for the three random walks, we see that a priori a site 
is colored green with probability about 59% in the first weighting, 30% with the 
symmetric weighting and only 7% in the reversed weighting. Yet, from a conditional 
point of view, the three colorings are identical. 

It turns out that the hinge is one of seven minimal graphs each of which can 
produce a phase transition on T r for any r > 2. The graphs are pictured in Plate 7. 
We say that a graph H is fertile if Hom(T r , H) has more than one simple, invariant 
Gibbs measure for some r and A; otherwise it is sterile. The fertile graphs are 
exactly those which contain one or more of the seven baby graphs in Plate 7 as an 
induced subgraph. It turns out that the value of r does not come into play: if the 
constraint graph is rich enough to produce a phase transition on any T r , then it 
does so for all r > 2. One way to state the result is as follows: 

Theorem 6.1 0] Fix r > 1 and let H be any constraint graph. Suppose that H 
satisfies the following two conditions: 

(a) Every looped node of H is adjacent to all other nodes of H ; 

(b) With its loops deleted, H is a complete multipartite graph. 

Then for every activity vector on H , there is a unique invariant Gibbs measure 
on the space Hom(T r ',iJ). 

If H fails either condition (a) or condition (b) then there is a set of activities 
A on H for which Hom(T r , H ) has at least two simple, invariant Gibbs measures, 
and therefore A can be obtained by more than one branching random walk. 

The proof of Theorcm l6.1l is far too complex to reproduce here, but reasonably 
straightforward in structure. First, a distinct pair of weightings yielding the same 
activity vector must be produced for each of the seven baby fertile graphs, and 
for each r > 2. Second, it must be demonstrated that if H contains one of the 
seven as an induced subgraph, then there are weightings (whose restrictions are 
close to those previously found) which induce phase transitions on Hom(T r , if). 
Third, a monotonicity argument is employed to show that if H satisfies conditions 
(a) and (b) of the theorem, then the map from w to A is injective. Finally, an easy 
graph-theoretical argument shows that H satisfies (a) and (b) precisely if it does 
not contain any of the seven baby fertile graphs as an induced subgraph. 
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7. An application to statistical physics 

Theorem 16.11 has many shortcomings, applying as it does only to hard con- 
straint models on the Bethe lattice, and we must also not forget that it considers 
only the very nicest Gibbs measures. The constraint graph of the hard-core model 
is sterile, yet it can have multiple Gibbs measures on T r (or, as we saw, on J?) if 
we relax the invariance condition. 

For H = the hinge, however, and for any r and A, there are multiple Gibbs 
measures on Hom(T r ,iJ) if and only if there are multiple simple, invariant Gibbs 
measures. Like the Ising model, the Widom-Rowlinson model exhibits spontaneous 
symmetry-breaking; indeed its relationship to the Ising model parallels the relation 
between nodes and edges of a graph. 

One of the nice properties known for the Ising model is that it can exhibit at 
most one critical point; but the proof of this fact does not work for the Widom- 
Rowlinson model. Indeed, in [2j the methods above are used to construct a board 
G for which Hom(G, H) has three (or more) calculable critical points, with H the 
hinge. Set A = A grcon = A rcc i, fixing A yc n ow = 1, so that the single parameter A 
controls the Widom-Rowlinson model. Then: 

Theorem 7.1 There exist < Ai < A2 < A3 and an infinite graph G, such that 
the Widom-Rowlinson model on G with activity A has a unique Gibbs measure for 
A G (0, Ai] U [A2, A3], and multiple Gibbs measures for A G (Ai, A2) U (A3, 00). 

The board G constructed in [2] is a tree, but not quite a regular one; it is made 
by dangling seven new pendant sites from each site of T . Readers are referred to 
that paper for the calculations, but the intuition is something like this. 

For low A the random coloring of G is mostly yellow, but as A rises, either 
green or red tends to take over the interior vertices as in T 40 . Then comes the third 
interval, where the septuplets of leaves, wanting to use both green and red, force 
more yellow on the interior vertices, relieving the pressure and restoring green-red 
symmetry. Finally the activity becomes so large that the random coloring is willing 
to give up red-green variety among the septuplets in order to avoid yellow interior 
vertices, and symmetry-breaking appears once again. 

It turns out that multiple critical points can be obtained for the hard-core 
model in a similar way. 

8. Dismantlable graphs 

In addition to the fertile and sterile graphs, a second graph dichotomy ap- 
pears repeatedly in our studies: dismantlable and non-dismantlable graphs. Co- 
incidentally, the term "dismantlable" as applied to graphs was coined by Richard 
Nowakowski and the second author |16j almost twenty years ago in another context 
entirely: a pursuit game on graphs. 

Two players, a cop C and a robber 1Z, compete on a fixed, finite, undirected 
graph H. We will assume that H is connected and has at least one edge, although 
the concepts make sense even without these assumptions. The cop begins by placing 
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herself at a node of her choice; the robber then does the same. Then the players 
alternate beginning with C, each moving to an adjacent node. The cop wins if she 
can "capture" the robber, that is, move onto the node occupied by the robber; 1Z 
wins by avoiding capture indefinitely. In doing so 1Z is free to move (or even place 
himself initially) onto the same node as the cop, although that would be unwise if 
the node were looped since then C could capture him at her next move. 

Evidently the robber can win on any loopless graph by placing himself at the 
same node as the cop and then shadowing her every move; among graphs in which 
every node is looped, C clearly wins on paths and loses on cycles of length 4 or 
more. (In the game as defined in |161 H"8] , there is in effect a loop at every node of 
H-) ^ 

The graph on which the game is played is said to be cop-win if C has a winning 
strategy, robber-win otherwise. The following structural characterization of cop-win 
graphs is proved in ^H] for the all-loops case, but in fact the proof (which is not 
difficult, and left here as an exercise) works fine in our more general context. 

Let N(i) be the neighborhood of node i in H and suppose there are nodes i 
and j in H such that N(i) C N(j). Then the map taking i to j, and every other 
node of H to itself, is a homomorphism from H to H\{i}. We call this a fold of the 
graph H . A finite graph H is dismantlable if there is a sequence of folds reducing 
H to a graph with one node (which will necessarily be looped). 

Note that dismantlable graphs are easily recognized in polynomial time. Plate 
8 shows some dismantlable and non-dismantlable graphs. 

The following theorem, from jjjj, collects a boatload of equivalent conditions. 

Theorem 8.1 The following are equivalent, for finite connected graphs H with at 
least one edge. 

1. H is dismantlable. 

2. H is cop-win. 

3. For every finite board G, Hom(G, H ) is connected. 

4- For every board G, and every pair tp, ip S Hom(G, H) agreeing on all but 
finitely many sites, there is a path in Hom(G, H) between tp and ip. 

5. There is some positive integer m such that, for every board G, every pair of 
sets U and V in G at distance at least m, and every pair of maps tp,ip £ 
Hom(G, -ff), there is a map 9 G Hom(G, H ) such that 9 agrees with tp on U 
and with ip on V . 

6. For every positive integer r, and every pair of maps tp, tp G Hom(T r , H ), there 
is a site u in T r with tp(u) ijj(u), a patch U containing u, and a map 
9 G Hom(T r , H) which agrees with ip on T r \ U and with ip on u. 

7. For every board G and activity vector X, if /1 is a measure on Hom(G, H ) 
satisfying the one-site condition, then (i is a Gibbs measure. 

8. For every finite board G and activity vector X, every stationary distribution 
for the point process V(G, H, A) is a Gibbs measure. 

9. For every board G of bounded degree such that Hom(G, H) is non-empty, 
there is an activity vector X such that there is a unique Gibbs measure on 
Hom(G,iT). 
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10. For every r, there is an activity vector A such that there is a unique Gibbs 
measure on Hom(T r ,lf). 

We will prove here what we think, to a graph theorist, is the most interesting of 
these equivalences — (i) and (Hi). Recall that two maps in Hom(G, H) are adjacent 
if they differ on one site of G. 

Let us first assume H is dismantlable. If it has only one node the connectivity 
of Hom(G, H) is trivial, since it has at most one element. Otherwise there are nodes 
i 7^ j in H with N(i) C N(j) and we may assume by induction that Hom(G, H') is 
connected for H' := H\ {i}. 

Define, for <p in Hom(G, H), the map <p' in Hom(G, H') (and also in Hom(G, H)) 
by changing all i's to j's in the image. If a and j3 are two maps in Hom(G, H) then 
there are paths from a to a' , a' to (3' and j3' to /3; so Hom(G, H) is connected as 
claimed. 

For the converse, let H be non-dismantlablc, and suppose that nonetheless 
Hom(G, H) is connected for all finite boards G; let q — \H\ be minimal with respect 
to these properties. 

If there are nodes i and j of H with N(i) C N(j), then H := H \ {i} is 
also non-dismantlable. In this case, we claim that the connectivity of Hom(G, H) 
implies connectivity of Hom(G, H'). To see this, define, for ip 6 Hom(G, H), the 
map (p' G Hom(G, H') by changing all i's to j's in the image as before. If a and 
[3 are two maps in Hom(G, H'), then we may connect them by a path (fx, . . . , ipt 
in Hom(G,-ff); now we observe that the not-necessarily distinct sequence of maps 
ip'x, . . . , tp' t connects a and f3 in Hom(G, H'). This contradicts the minimality of i?, 
so we may assume from now on that there is no pair of nodes i ^ j in H with 
N(i) C N(J). 

Now let G be the 'weak' square of H, that is, the graph whose nodes are 
ordered pairs (ii, ia) of nodes of H with (11,12) ~ (ji, 32) just when i\ ~ ji and 
*2 ~ J2- There are two natural homomorphisms from G to -ff, the projections tti and 
7T2, where 7Ti(ii, i%) = i% and ^(ii, 12) = «2! we claim that 7Ti is an isolated point of 
the graph Hom(G,iJ), which certainly implies that Hom(G,7J) is disconnected. 

If not, there is a map it' taking (say) 12) to k ^ i\ and otherwise agreeing 
with 7i"i . Let J2 be a fixed neighbor of «2 and ji any neighbor of i\. Then {i\, 12) ~ 
(ill J2), and hence fc ~ ji in i?. We have shown that every neighbor of i\ is also a 
neighbor of fc, contradicting the assumption that no such pair of nodes exists in H. 
This completes the proof. 

For the last part of the proof, there is also a simpler (and smaller) construction 
that works provided H has at least one loop: see [5] or Cooper, Dyer and Frieze [5]. 

We have seen that, for a dismantlable constraint graph H, and any board G of 
bounded degree, there is some A (which can be taken to depend only on H and the 
maximum degree of G) such that there is a unique Gibbs measure on Hom(G, H). 
Dyer, Jerrum and Vigoda |10| have proved a "rapid mixing" counterpart to this 
result: given a dismantlable H, and a degree bound A, there is some A such that 
the point process V(G, H, A) is rapidly mixing for all finite graphs G with maximum 
degree at most A. Of course, if H is not dismantlable, then no such result can be 
true as Hom(G, H) need not be connected. 
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9. Random colorings of the cayley tree 

We have observed that ordinary "proper" g-colorings of a graph G arc maps 
in Hom(G, K q ); since K q is sterile, there is never more than one Gibbs measure for 
G = T r . However, we see even in the case q = 2 that multiple Gibbs measures 
exist, because each of the two 2-colorings of T r determines by itself a trivial Gibbs 
measure, as does any convex combination. However, only the \, \ combination is 
invariant under parity-changing automorphisms of T r . 

All of the Gibbs measures in the q = 2 case are, however, simple and invariant 
under all the parity-preserving automorphisms of the board. Such Gibbs measures 
are neededed to realize the phase transition for the hard-core model as well, so it 
is not surprising that it is useful to relax our requirements slightly and to consider 
these semi-invariant simple Gibbs measures. 

Fortunately we don't have to throw away all our work on invariant Gibbs 
measures in moving to semi-invariant ones. Given a constraint graph H on nodes 
1,2,..., q which is connected and not bipartite, we form its bipartite "double", 
denoted 2H, as follows: the nodes of 2H are {1,2..., q}U {— 1, —2, . . . , — q} with an 
edge between i and j just when i <~ —j or — i <~ j in H. Note that 2H is loopless; 
a loop at node i in H becomes the edge {—i, i} in 2H. 

A homomorphism i\> from T r to 2H induces a homomorphism \ip\ to H via 
\i>\(v) — \ip{v)\- In the reverse direction, a map <p in Hom(T r ,i/) may be trans- 
formed to a map (p in Hom(T r , 2H), by putting <p{v) — ip(v) for even sites v e T r 
and <p(v) = —ip(v) for odd v. 

Let A = (Ai, . . . , A„) be an activity vector for H and suppose that fi is a simple 
invariant Gibbs measure on Hom(T r , H) corresponding to A. From [i we can obtain 
a simple invariant Gibbs measure ft on Hom(T r ,2£f) by selecting ip from ji, and 
flipping a fair coin to decide between (p (as defined above) and —(p. Obviously p, 
yields the activity vector A on 2H given by Ai = A|j|. Furthermore, the weights on 
H which produce /i extend to 2H by w-i = Wi. 

Conversely, suppose v is a simple invariant Gibbs measure on Hom(T r , 2i?) 
whose activity vector satisfies A_i = Aj for each i. Then the measure \v\, obtained 
by choosing ip from v and taking its absolute value, is certainly an invariant Gibbs 
measure on Hom(T r , H) for A |"{1, . . . , q}, but is it simple? 

In fact, if the weights on 2H which produce v do not satisfy w^i — cwi, then 
\v\ will fail to be simple. To sec this, observe that if the weights are not proportional 
then there are nodes i <~ j of H such that P-i,-j ^ Pij in the random walk on 2H. 
Suppose that l^l is conditioned on the color of the root w of T r being fixed at i, 
and let x and y be distinct neighbors of w. Set a = Vi{ip{w) = i \ |^>(tu)| = i). 
Then 

Pi(\ip\(x) = j) = (1 - a)p-i-j + apij 

but 

Pi(Mx) = j A |V|(y) = j) = (1 - tt> 2 ,,-j + a Ph > P <M X ) = J) 2 
so the colors of x and y are not independent given 
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However, we can recover simplicity at the expense of one bit worth of symme- 
try. Let v + be v conditioned on tp(u) > 0, and define v~ similarly. Then \v Jr \ and 
\v~\ are essentially the same as v + and v~ , respectively, and all are simple; but 
these measures are only semi-invariant. 

On the other hand, suppose fi is a simple semi-invariant Gibbs measure on 
Hom(T r , H). Let be a parity-reversing automorphism of T r and define // := fj,o9, 
so that \[i + \\J is fully invariant (but generally no longer simple). However, 
v := |/2+ |(— //) is a simple and invariant Gibbs measure on Hom(T r ,2_ff), thus 
given by a node- weighted random walk on 2H. We can recover /x as v +1 hence: 

Theorem 9.1 Every simple semi- invariant Gibbs measure on Hom(T r ,H) is ob- 
tainable from a node-weighted branching random walk on 2H , with its initial state 
drawn from the stationary distribution on positive nodes of 2H . 

Suppose, instead of beginning with a measure, we start by weighting the nodes 
of 2H and creating a Gibbs measure as in Theorem 19. II Suppose the activities of 
the measure are {Xi : i = ±1, . . . ,±g}. By identifying color —i with i for each 
i > 0, we create a measure on H-colorings, but this will not be a Gibbs measure 
unless it happens that (A_i, . . . , X- q ) is proportional to (Ai, . . . , X q ). 

We could assure this easily enough by making the weights proportional as 
well, e.g. by W-i — Wi\ then the resulting measure on Hom(G,-ff) could have been 
obtained directly by applying these weights to H, and is thus a fully invariant 
simple Gibbs measure. To get new, semi-invariant Gibbs measures on Hom(G, H ) , 
we must somehow devise weights for 2H such that u>-i Wi yet X-i oc A^. 

Restated with slightly different notation, simple semi-invariant Gibbs measures 
are in 1-1 correspondence with solutions to the "fundamental equations" 



for i = 1, . . . , q. Such a solution will be invariant if Ui — Vi for each i. 

Plate 9 illustrates a semi-invariant, but not invariant, simple Gibbs measure 
for uniform 3-colorings of T 3 . Approximate weights of the nodes of 2H = 2K% are 
given along with part of a sample coloring drawn from this measure. Additional 
measures may be obtained by permuting the colors or by making all the weights 
equal (invariant case). 

Results for g-colorings of TP', with q > 2 and r > 1, are as follows: 

When q < r + 1 , all choices of activity vector including the uniform case yield 
multiple simple semi-invariant Gibbs measures. 

When q > r + 1 , there is only one simple semi-invariant Gibbs measure for 
the uniform activity vector, but multiple simple semi-invariant Gibbs measures for 
some other choices of activity vector. 

The critical case is at q = r + 1, that is, when the number of colors is equal 
to the degree of the Cayley tree. Here it turns out that there are multiple simple 
semi-invariant Gibbs measures for all activity vectors except the uniform case, where 
there is just one. 
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When q > r + l and the activities are equal, the unique simple semi- invariant 
Gibbs measure is in fact the only Gibbs measure of any kind. This was conjectured 
in |S] but proved only for q > cr, with fixed c > 1; Jonasson ^3] nas recently, 
and very nicely, finished the job. Jonasson's result is in a sense a special case of 
the conjecture that the Markov chain of q-colorings of a finite graph of maximum 
degree less than q, which progresses by choosing and recoloring sites randomly one 
at a time, mixes rapidly. So far the best result is Vigoda's [2U] which proves this if 
the maximum degree is at most 6g/ll. 

When q < r+l there are lots of other Gibbs measures, including ones we call 
frozen. These come about because it is possible for a measure to satisfy the Gibbs 
condition in a trivial and somewhat unsatisfactory way. For example, suppose we 
are ^-coloring TP' (with root w) for some q < r+l, and let tp be any fixed coloring in 
which the children of every node exhibit all colors other than the color of the parent. 
Let fi be the measure which assigns probability 1 to ip. Then for any finite patch 
U, which we can assume to be a subtree including the root, the colors on dU force 
the colors on the leaves of U, and we can continue inwards to show that the original 
coloring ip \ U is the only one consistent with the colors on dU. Thus \i satisfies 
the Gibbs condition trivially, and is also vacuously simple — but not invariant or 
semi-invariant. We call a Gibbs measure of this type "frozen". A frozen state of 
Hom(T 2 , i^) is illustrated in Plate 10. For more about frozen Gibbs measures the 
reader is referred to 0. 

In a soft constraint model such as the Potts model, frozen Gibbs measures can 
only occur at zero temperature. Since most of the time statistical physicists are in- 
terested only in phases which exist at some positive temperature (and have positive 
entropy), frozen measures are generally absent from the statistical physics litera- 
ture. However, they are interesting combinatorially and motivate some definitions 
in the next section. 

10. Prom statistical physics back to graph theory 

We conclude these notes with a theorem and a conjecture in "pure" graph 
theory, stripped of probability and physics, but suggested by the many ideas which 
have appeared in earlier sections. 

Suppose H is bipartite and we are given some ip G Hom(T 1 , H) , where the 
sites of T 1 are labeled by the integers Z. Then knowing ip(n) even for a very large 
n tells us something about ip(0), namely which "part" of H it is in. We call this 
phenomenon long range action, and define it on Cayley trees as follows: If there 
is a ip £ Hom(T r , H) and a node i £ H such that for any n, no ip G Hom(T r , H) 
agreeing with tp on the sites at distance n from the root can have spin i at the root, 
we say Hom(T r , H) has long range action. 

Theorem 10.1 If H is k-colorable then Hom(T fe_1 , H) has long range action. 

For example, the coloring described at the end of the previous section, which 
gives rise to a frozen Gibbs measure, shows that Hom(T 2 ,ii3) (more generally, 
Hom(T r , K r+ i)) has long range action. We also see from Theorem 18.11 that 
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Hom(T r ,ff) has long range action for no r if and only if H is dismantlable; of 
course then H has at least one loop and therefore has infinite chromatic number. 

Note that the theorem connects a statement about homomorphisms from H 
to a statement about homomorphisms to H . However, it is difficult to see how to 
turn a fc-coloring of H into a suitable map in Hom(T fc_1 , H). Suppose, for instance, 
that H is the 5-cycle C5, with nodes represented by the integers modulo 5. We can 
get a completely frozen map in Hom(T 2 ,C5) by making sure we use both i + 1 and 
i— 1 on the children of any site of spin i. But what has this map got to do with any 
3-coloring of C5? 

The proof of Theorem llO.il found in |7j, uses a vector- valued generalization 
of coloring to construct the required map in Hom(T fc_1 , H). 

We now move from long range action to the familiar notion of connectivity. 
Theorem l8.ll — in fact, the part whose proof is given above — tells us that Hom(G, H) 
is connected for any finite G just when H is dismantlable. Suppose we restrict 
ourselves to boards of bounded degree? If, for example, H is bipartite, Hom(K2, H) 
is already disconnected. If H = K q then Hom(K q , H) is extremely disconnected, 
consisting of nl isolated maps. By analogy with Theorem llO.il we should perhaps 
be able to prove: 

Conjecture 10.2 If H is k-colorable then Hom(G, H) is disconnected for some 
finite G of maximum degree less than k. 

A proof for k — 3 appears in [7] and Lovasz ^1] has shown that the conjecture 
holds for k = 4 as well. We think that a proof of Conjecture 110.21 would have to 
capture some basic truths about graphs and combinatorial topology, and fervently 
hope that some reader of these notes will take up the challenge. 
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